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Abstract 

We investigate the radiative dileptonic decays B s (Bd) — > within the standard 

model. Using the constituent quark model, the branching ratios turn out to be around 
5 x 10~ 9 for B s — ► 7/i + /x~ and around 6 x 10~ 10 for — ► 7/U + /i~, with slightly larger 
values for B s (Bd) — > 7e + e _ . The differential rate as a function of the dilepton invariant 
mass is given. The possibility of using these processes to determine the decay constants 
of B s or B<i is discussed. 
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1 Introduction 



Rare decays induced by flavor changing neutral currents can be used as tests of the Standard 
Model (SM), and are sensitive to new physics. It is expected that at future B factories and 
fixed target machines, other rare decay channels of the bottom quark |1J will be discovered in 
addition to the observed b — > S7 transition 0. These processes also offer useful information 
for extracting the fundamental parameters of the SM, such as \V u b\ 0. 

Rare decays can also serve as alternative channels to measure some elementary hadronic 
parameters. For instance, the decay constants fs q , q = s,d can be extracted from B q — > juP 
[Q. In the present work, we investigate the possibility of using the channel, B q — ► 
to determine the decay constants. Pure leptonic decays of heavy pseudoscalar mesons into 
light lepton pairs are helicity suppressed, their branching ratios are of the order of 10~ 9 for 
B s — > , and 10 -14 for B s — > e + e~ ||. This makes it difficult to determine fs s from 
these processes. For Bd the situation is even worse due to the smaller CKM angle. If a 
photon is emitted in addition to the lepton pair, then the mechanism of helicity suppression 
will not hold any longer and larger branching ratios are expected. Although the process 
B s — > t + t~ is free from helicity suppression and its branching ratio is around 8 x 10 -7 in the 
SM , it is likely to be compatible with the decays into lepton pairs only when its efficiency 
is better than 10~ 2 . 

In Section 2 the relevant effective Hamiltonian will be given in the SM, and the constituent 
quark model will be used to give the numerical predictions. Finally, we make some comments 
in Section 3. 

2 Model calculations 

The most important contributions to B q — > (/ = e, /i) stem from the effective Hamil- 

tonian which induces the pure leptonic processes B g — > l + l~ . Let us start with the short 
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distance contributions to the latter one, which include box, Z and photon penguin diagrams. 
The Feynman diagrams are displayed in Fig.l. The QCD corrected effective Hamiltonian in 
the SM is % g: 



aG 



Heff ~ 7^ VtbK 



C e 9 ff qYP L b + fy»P R b\ h»l + C 10 (q^P L b)l llxl5 l 



with P L — (1 — 75)/2, P R — (1 + 7s)/2, and p = p + + p_ is the momentum of the lepton 
pair. The QCD corrected Wilson coefficients CV(m;,), Cg^(rrib), and Cio(mb) are given by 
Misiak §. 

Because of the lightness of the leptons e and /x, the pure leptonic processes B q — » Z + Z~ are 
suppressed by helicity. If a photon line is attached to any of the charged lines in Fig.l, the 
situation will be different: helicity suppression is overcome. When the photon is attached to 
internal lines, there will be a suppression factor of m^/M^r in the Wilson coefficient, since 
the resulting operators(dimension-8 ) are two orders higher in dimensions than the usual 
ones(dimension-6) . 

One may expect that the additional two diagrams displayed in Fig. 2 (a) and (b) might 
also contribute to the decay B q — > However, our calculation shows that the contri- 

bution of Fig. 2(b) is quite small and its influence can be neglected in the numerical results. 
Fig. 2 (a) involves a singularity, since the intermediate s quark can be on its mass shell. The 
influence of this singularity is an artifact of the constituent quark model employed here. If 
we calculate the contribution of Fig. 2(a) to the differential decay rate anywhere except near 
the singularity we find it to be negligible ( at most a few percent). Therefore we neglect it 
altogether. 

There are also long distance contributions to B q — > l + l^'~f. For instance, in the B s decay, 
there are cascade processes at the hadronic level: 

b s - 0(0', JM-h, 0(0', J/V-) - i + r- (2) 



These involve the on-shell 0, <p', J/ip... and are similar to b — > sJ/ip — > sl + l They 



are estimated to be of order of 10 9 , which is about the same order of the short distance 



contribution. All these long distance effects, which are shown as sharp peaks in the invariant 
mass spectrum of the lepton pair, will not be included in the following. 

Now using a simple constituent quark model (see, for example, ||), we calculate the 
amplitude for diagrams with photon emitted from external fermion lines. First, for diagrams 
with photon attached to external lepton lines, we get the amplitude proportional to: 

K+ fa - m i 



qJv.b L e v l 
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(P- +Pj) 2 ~ m i 
Using the decay constant definition: 



(p+ + P7) 2 - m « 



'7 



<0\q7^ 5 b\B>=-f Bq p B , 



I. 

(3) 



(4) 



we can easily calculate the results of eqn.(|3|). Note that, there is no contribution from the 
operator O7 in this circumstance, since the definition of the decay constant shows that: 



< 0\qa^P R b\B >= 0. 



(5) 



For the other two diagrams with photon emitted from the external quark lines (b or q), 
the amplitude is proportional to: 



(Pq ■ Pi) 



(Pb ■ P 7 ) 



C e J f l ltl l + c w i w i) 



+ 2 — — e v q 



k~ fa + m b v v fa- fa + m q 



-7+7 



b («TmO- 



(6) 



(Pb'Pj) (Pq-Pj) 

In the constituent quark model, p^ = (m q /m B )p B , and p% = {m h /m B )p B ^. Then applying 
eqn.(H), neglecting terms suppressed by m q /rrib, q — d,s, we get: 



A — C it a i3^ v e c ^p 1 p l Bq + (p 7At e 7!/ — p lv e lil )p Bq 
with 



' c eff _ 2C 7 m2 B q 



P 



(7) 



c = 



eaG F f B m B 



-v tb v; q . 



l2\/2iim q {j) Bq -p 7 ) 

After squaring the total decay amplitude, we find that the contribution of diagrams with 
photon attached to the lepton lines (eqn.(|) ) is proportional to mi, which is negligible in 



numerical calculations. The main contribution to B q — > 7/+/" is thus from the square of 
eqn.(0). Performing the phase space integration over one of the two Dalitz variables, we get 
the differential decay width versus s = (p + + p„) 2 / m Bq : 

We see that the differential rate is proportional to (fB q / m q) 2 - Using a = 1/137, m s = 0.51 
GeV, m t = 176 GeV and = 0.04, the differential decay rate as a function of s is 

displayed in Fig.3, which shows that the contribution from soft photons, corresponding to 
large s region, is negligibly small. 



After integration, we also get: 
T(B S -> 7 e + e- 

T(B 8 -> 



3.0 x 10 



2.2 x 10 



-21 



-21 



X 



0.2GeV, 
0.2GeV 



GeV. 



GeV. 



(9) 



For B d meson decay, we take |V#,V^| = 0.01 and m d = 0.35 GeV. The decay widths are then: 

\ 2 



T{B d 

r(B d - 



7e e 



3.6 x 10 



2.7 x 10 



-22 



-22 



B,, 



0.2GeV / 



If the lifetimes are taken as t(B s 



\0.2GeV / 
1.34 x 10~ 12 s, r{B d ) = 



GeV, 



GeV. 



(10) 



1.50 x 10" 12 s |g, and if 



/s = 200MeV is used, they correspond to the branching ratios: 



B(B S - 


-> e + e 7) 


= 6.2 x 10~ 9 , 




B{B a - 




= 4.6 x 10~ 9 , 
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B(B d - 


-> e + e~7) 


= 8.2 x 10~ 10 , 




B(B d - 




= 6.2 x 10~ 10 . 





Note that the branching ratios we get for radiative decays are just the same order of the 

pure leptonic decay B q — ► . The decay rates for pure leptonic decays are: 

a 2 G 2 F flm Bq mf 



T (B q - rr) 



16vr 3 



v tb v t * q \ 2 c 2 



10 • 



(12) 



Using the same parameters as in the radiative decay, the numerical results are: 

B(B S -> e + e~) = 6.1 x 1(T 14 , 

B(B S -> //V") = 2.6 x KT 9 , (13) 
B(B d -> e + e~) = 4.2 x 1(T 15 , 
5(5^ -> = 1-8 x 1(T 10 . 

So for the e + e~ channel, the radiative decay dominates over the dileptonic one, while for the 
fj, + fj,~ channel, the branching ratio of the radiative decay is a little larger than that for the 
pure leptonic decay. 

3 Conclusions 

We predict the branching ratios in the SM for B s — > ^yl + l~ to be around 5 x 1CT 9 for 
B s — > 7/i + /i _ and around 6 x 10~ 10 for Ba — > 7/i + /i~, with slightly larger values for B s (Bd) — > 
7e + e~. With these predictions, they will hopefully be detected at LHC-B. thus provide 
alternative channels for measuring fB a / m s (/s d /^d)- In LHC-B, approximately 6 x 10 11 
(2 x 10 11 ) Bd (B s ) mesons are expected per year, therefore there is a good chance of observing 
the decays considered here. 

We conclude with comments on the experimental extractions of the decay constants. In 
the measurements of the pure leptonic decays, processes with additional soft photons emitted 
(say, photons with energy less than 50 MeV) might contribute due to the inability of the 
detector to separate them. From our calculation, we find that these soft photon processes are 
not important in the radiative dileptonic decays, see Fig.3. Since the efficiency of detecting 
photons may be low, it will be useful for experiments to have the sum of the radiative decay 
and the pure leptonic processes: 

B(B S -> /!+/!- + 7/i + /i") = 7.3 x 1(T 9 , 

B(B d -> yuV" + 7At + /0 = 8.1 x 1(T 10 . (14) 
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The total n + fi (7) branching ratios for B s and decay are almost the same as that of the 
e + e~7 case. 
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Figure Captions 

Fig.l Feynman diagrams in standard model 

Fig. 2 Additional diagrams which contribute 
ators are O7 (black dot). 

Fig. 3 Differential decay rates of B s — > 



for bq — > l + l . 

to B s — > / ~fl + l~ . The relevant penguin oper- 
versus s = (p + + p-) 2 /m 2 Bg . 
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